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Abstract 

The line of phase transitions, separating the confinement and the 
q_) \ Coulomb phases in the four-dimensional pure compact U(l) gauge 

theory with extended Wilson action, is reconsidered. We present new 
numerical evidence that a part of this line, including the original Wil- 
^ | son action, is of second order. By means of a high precision simu- 

;_, ' lation on homogeneous lattices on a sphere we find that along this 

line the scaling behavior is determined by one fixed point with dis- 
tinctly non-Gaussian critical exponent v = 0.365(8). This makes the 
existence of a nontrivial and nonasymptotically free four-dimensional 
pure U(l) gauge theory in the continuum very probable. The univer- 
sality and duality arguments suggest that this conclusion holds also 
for the monopole loop gas, for the noncompact abelian Higgs model at 
large negative squared bare mass, and for the corresponding effective 
string theory. 



1 Introduction 



In contradistinction to lower dimensions, the only firmly established quantum 
field theories in four dimensions (4D) are either asymptotically free or so- 
called trivial theories. Both are defined in the vicinity of Gaussian fixed 
points, i.e. of noninteracting limit cases. In spite of numerous suggestions 
and circumstantial evidence, until now no candidate for a non-Gaussian fixed 
point in 4D has been established. 

The suitability of numerical simulations on the lattice for a confirmation 
of the existence of non-Gaussian fixed points, and for an investigation of 
their properties, has been demonstrated in dimensions lower than four in nu- 
merous applications. For example, non-Gaussian values of critical exponents 
can be determined by means of finite-size scaling (FSS) or renormalization 
group (RG) analysis. Several attempts to use a similar approach in 4D have 
encountered severe problems, however. 

In this letter we reconsider the oldest candidate for a non-Gaussian fixed 
point in the 4D lattice field theory, the phase transition between the con- 
finement and the Coulomb phases in the pure compact U(l) gauge theory 
PU |2] with Wilson action and extended Wilson action. After various pio- 
neering studies, e.g. [3IU|3], the more detailed investigations were hindered 
mainly by a weak two-state signal (HI Ej. This obscures the order of the 
phase transition and makes it difficult to determine critical exponents reli- 
ably. We demonstrate that the problems encountered, when considering the 
continuum limit at this phase transition, can be surmounted. The clues are 
the observation that the two-state signal disappears on lattices with sphere- 
like topology, the construction of homogeneous spherical lattices, the use of 
modern FSS analysis techniques, and larger computer resources. 

We find that at the confinement - Coulomb phase transition at strong bare 
gauge coupling, g = 0(1), the model exhibits a 2 nd order scaling behavior 
well described by the values of the correlation length critical exponent v in 
the range v = 0.35 — 0.40. The measurements have been performed at various 
couplings and by different methods. The most reliable determination gives 

v = 0.365(8). (1) 

These results are quite different from v = 0.25, expected at a 1 st order 
transition, as well as from v = 0.5, obtained in a Gaussian theory or in 
the mean field approximation. This strongly suggests the existence of a 
continuum pure U(l) gauge theory with properties different from theories 
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governed by Gaussian fixed points with or without logarithmic corrections. 
It can be obtained from the lattice theory by the RG techniques. 

Detailed numerical evidence for these claims will be presented elsewhere |S] . 
Some preliminary results have been published in Refs. [§1 110j. 

Since the pure U(l) lattice gauge theory with the Villain (periodic Gaus- 
sian) action presumably belongs to the same universality class [TT], rigorous 
dual relationships imply that also the following 4D models possess a contin- 
uum limit described by the same fixed point: the Coulomb gas of monopole 
loops[T2], the noncompact U(l) Higgs model at large negative squared bare 
mass (frozen 4D superconductor) |T3*1 IT3j . and an effective string theory 
equivalent to this Higgs model[T5]. 

These findings raise once again the question, whether in strongly inter- 
acting 4D gauge field theories further non-Gaussian fixed points exist, that 
might be of use for theories beyond the standard model. 

2 Earlier results 

The pure compact U(l) gauge theory on the 4D cubic lattice with periodic 
boundary conditions (4D torus) can be described by the extended Wilson 
action [S] 

S = -^w P [/3cos(e P ) + 7 cos(2e P )] , (2) 
p 

with wp = 1. Here 0p G [0, 2ti) is the plaquette angle, i.e. the argument of 
the product of U(l) link variables along a plaquette P. Taking 6p = a 2 gF ilV) 
where a is the lattice spacing, and (3 + 4 7 = l/g 2 , one obtains for weak 
coupling g the usual continuum action S = \ J d A xF 2 u . 

In one of the very first studies ||, restricted to the 7 = case (Wilson 
action pP) and small lattices, a behavior consistent with a 2 nd order phase 
transition at (5 ~ 1 was observed, and v ~ 1/3 was found. However, any 
inference about the continuum limit has been hindered by the subsequent 
observation of a two-state signal on larger, but finite lattices^. This could 
either be a finite-size effect, or it could imply that the phase transition at 
7 = is actually of 1 st order, preventing a continuum limit at 7 = 0. 

In the model with the extended Wilson action (J2J), it was found that 
the confinement-Coulomb phase transition is clearly of 1 st order for 7 > 
0.2, and weakens with decreasing 7 EI- Various studies suggested that 
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the transition becomes 2 nd order at slightly negative 7 0, or around 7 = 

opiiin!. 

The order of the transition at 7 = has remained a controversial subject 
[T%1 Uni 1201 EE] Though the values v ~ 0.3 — 0.4 have been obtained consis- 
tently by various methods jU HH1 H3 122 I2H] , the continuum limit has not 
appeared to be possible. 

Also the hope that the continuum limit might be taken at least at negative 
7 was spoiled by the observation of a weaker, but still significant, two-state 
signal on finite lattices even there |7j. Though this signal is probably only a 
finite-size effect, and the transition in the infinite volume limit is genuinely 
of 2 nd order, it impedes a precise FSS and RG analysis. 

It was known that monopole loops winding around the toroidal lattice oc- 
cur |2*4*] ITo] and cause difficulties in simulations with local update algorithms. 
Suspecting that this might be a reason for the two-state signal, two of the 
present authors performed simulations at 7 = 0, using the 4D surface of a 5D 
cubic lattice instead of the torus. They observed that on such lattices with 
spherelike topology the two-state signal vanishes [IE] . This suggests that the 
two-state signal at 7 < is related to the nontrivial topology of the toroidal 
lattice. 

Related observations have been made for the Schwinger model [25J. On 
the other hand, it has been checked in spin models, that weak two-state 
signals are not washed out on lattices with spherelike topology, if they are 
due to a genuine 1 st order transition [10J. 

However, the lattice on the surface of a cube is rather inhomogeneous and 
causes complex finite-size effects, preventing a reliable FSS analysis. 

3 Present methods 

For our present study at 7 < 0, we have chosen again lattices with spherelike 
topology. To alleviate the problem of inhomogeneity, we have used lattices 
obtained by projecting the 4D surface SH[N] of a 5D cubic lattice N 5 onto 
a concentric 4D sphere. On such a spherical lattice S[N], the curvatures 
concentrated on the corners, edges, etc., of the original lattice SH[N] are 
approximately homogenized over the whole sphere by the weight factors 

w P = A' P /Ap (3) 
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in the action (J2J). Ap and A' P are the areas of the plaquette P on S[N], and 
of its dual, as on any irregular, e.g. random lattice [2H] • These areas are 
determined by means of a two-triangle approximation of the plaquettes. 

It has been checked in some spin and gauge models with 2 nd order transi- 
tions that universality for spherical lattices holds, and that the FSS analysis 
works very well if V X I D is used as a linear size parameter, V = |X)pU>p 
being the volume of the sphere S[N] [TTH l2"7j. 

Another new feature in our study of the theory is the FSS analysis of the 
first zero z Q of the partition function in the complex plane of the coupling (3 
(Fisher zero). Applying the multihistogram reweighting method [2E] for its 
determination, the expected FSS behavior 

lmz Q <xV- 1,Du (4) 

has been used for measuring v. This has turned out to be superior to - 
though consistent with — the more common FSS analysis of specific-heat 
and cumulant extrema. 

Finally, we have performed the FSS study of the confinement- Coulomb 
transition not only at 7 = 0, but also at 7 = —0.2 and —0.5. This allows 
us to compare the behavior of the system at 7 = 0, where the order of the 
transition is disputed, with the commonly expected 2 nd order behavior at 
negative 7, and to test the universality of the critical properties. 

4 Results 

We have performed simulations [Hj on S[N] for N between 4 and 12. Note 
that S [12] has about (19. 6) 4 lattice points. The values of (3 were chosen 
in the immediate neighborhood of its critical values for 7 = 0, —0.2, and 
—0.5. For each lattice size at each 7 value, we have accumulated typically 
10 6 updates distributed over 8-12/3 points. 

We have found no indication of a two-state signal, neither in the individ- 
ual, nor in the multihistogram distributions of e = (J2p w p cos (@p))/(Sp Wp), 
at any of the three 7 values. This is demonstrated in Fig. 1 for 7 = 0. The 
values of the studied cumulants at their respective extrema are compatible 
with a 2 nd order transition. The critical behavior at all three 7 values is very 
similar, except that the transition weakens with decreasing 7, which means 
that larger lattices are needed for the same height of the specific-heat peak. 
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Figure 1: Histogram of the energy density e on the S [10] lattice (about 16 4 
points) at 7 = and (3 = 1.0180, very close to the pseudocritical point 
(3 pc = 1.01835(4) on that lattice, h is the relative occurence of e-values. 



In test runs at 7 = +0.2, a two-state signal has been clearly observed. 
This confirms that at sufficiently large positive 7 the phase transition is of 
1 st order, and that the spherical lattice S[N] does not wash out such a signal. 

Furthermore, at all three investigated values 7 < 0, the FSS analysis 
assuming a 2 nd order transition works remarkably well and leads to consis- 
tent results for all observables [5]. All our evidence thus points towards the 
conclusion that the phase transition is of 2 nd order for 7 < 0. 

In Table 1 we present results for the critical exponent v obtained from 
all the data with N > 4. The most reliable ones come from the FSS analysis 
of the Fisher zero (first column). The approximate agreement between the 
obtained values of v at all three 7 values demonstrates that the confinement- 
Coulomb critical line belongs to one universality class. The FSS behavior 
of Im^o according to (J3J), and the consistency of this behavior at different 
7 are illustrated in Fig. 2. The value (JI} has been obtained from a joint fit 
by means of (JIJ) to the data for Im^o with N > 6 (parallel straight lines in 
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Table 1: Results for v from fits to Imz , cy, n 
The indicated errors are statistical. 



x and V CL B,min at various 7. 



7 


Im z 


CV,max 


VcLB,min 


0. 


0.345(3) 


0.361(6) 


0.361(6) 


-0.2 


0.378(7) 


0.374(6) 


0.365(6) 


-0.5 


0.368(8) 


0.404(9) 


0.396(9) 



Fig. 2). 

The next two columns in Table 1 show the results from the FSS anal- 
ysis of the maximum of the specific-heat Cy >max and the minimum of the 
Challa-Landau-Binder cumulant VcLB,min- Here certain assumptions about 
nonleading terms have been necessary [Hj, and small systematic errors are 
therefore possible. Nevertheless, all the shown results are consistent with v 
lying in the interval 0.35 — 0.40, thus supporting universality. 
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Figure 2: Joint fit to Imzo according to (jlj) for all three 7 values and N > 6. 
The values of L = V 1 ^ correspond to iV = 4 — 10, 12. 
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5 Discussion 



The physical content of the continuum limit of the pure compact U(l) gauge 
theory at the confinement- Coulomb phase transition, discussed e.g. in |13| 
l2*Hj . depends on the phase from which the critical line is approached. In 
the confinement phase, a confining theory with monopole condensate is ex- 
pected, as the string tension scales with a critical exponent consistent with 
the value (JIJ[22|. The physical spectrum consists of various gauge balls. In 
the Coulomb phase, massless photon and massive magnetic monopoles, both 
already observed in Monte Carlo simulations [201 El], should be present. The 
renormalized electric charge e r is large but finite|32|l22|. and has presumably 
a universal value |H2J EUS E2j . The numerical result e 2 r /An = 0.20(2) [221 a S rees 
with the Liischer bound (Slj, as explained in [3T]- 

To our knowledge, the existence of such continuum quantum field the- 
ories in 4D is in no way indicated by the perturbation expansion. The 
non-Gaussian character of the fixed point might be rather understood as 
a consequence of the complex dynamics of the systems obtained in the dual 
representation of the theory with Villain action [T21 HB1 EH H51I3T] . 

There are some questions deserving further discussion. For 7 < 0, the 
studied theory does not satisfy reflection positivity, which is a sufficient, 
albeit not necessary condition for unitarity. If the phase transition in the 
reflection positive case at 7 = is of 2 nd order, as strongly suggested by our 
results, then unitarity should hold also for 7 < by universality arguments. 
If it is of weak 1 st order, unitarity at 7 < is made plausible by our finding 
that the scaling behavior at 7 = (on lattices of limited size), and at 7 < 
is the same, and that the regions with 7 < and 7 > are connected by the 
RG flows Pl|. 

Though improbable, it is not completely excluded that the phase tran- 
sition we have studied is of 1 st order for any 7 < 35J, but so weak, that 
the two-state signal cannot be detected by the currently available numer- 
ical means. Then the correlation length would remain finite for all finite 
7. Presumably, the effective cutoff, arising in this way, may be made ar- 
bitrarily large by suitably decreasing 7. However, unitarity then might be 
questionable. 

At small 7 > 0, where the order of the transition probably changes, a 
tricritical point with special values of indices is expected^. As the domain 
of dominance of such a point is unknown, it could be that the measured value 
of v, and the corresponding non-Gaussian fixed point are actually tricritical. 
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Alternatively, the change of the order might be more complicated than in 
metamagnets with tricritical points, since in our case only nonlocal order 
parameters are available. 

Though this may seem plausible, we cannot firmly conclude that the 
two-state signal observed on toroidal lattices is due to the monopole loops 
winding around them. The evidence for such an explanation is controversial 
(23J EH EI3 EJj , and the question, which configurations make the difference 
between the compact U(l) gauge theory on finite toroidal and spherical lat- 
tices, requires further study. The thermodynamic limit ought to be the same. 

We thank D.P. Landau, T. Lippert, M.I. Polikarpov, and U.-J. Wiese for 
discussions. The computations have been performed on the CRAY-YMP of 
HLRZ Jiilich and at the Parallel Compute Server of KFU Graz. J.J. was 
supported by Deutsches BMBF. 
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